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ABSTRACT 

The " density-density" correlation function of conduction electrons in metal 
is investigated. It is shown, that the asymptotic behaviour of the CF depends 
on the shape and the local geometry of the Fermi surface. In particular, the 
exponent of power law which describes the damping of Friedel oscillations at 
r — > oo(— 4 for an isotropic Fermi gas) is determined by local geometry of the 
FS. The applications of the obtained results to calculations of the CF in a metal 
near the electron topological transition and of the RKKY exchange integral are 
considered as well. 

Classification numbers: 71.45G, 71.25P. 
1. Introduction. 

In this paper we investigate the "density-density" correlation function (CF) 
of conduction electrons in a metal at T = . 

It is well-known, that correlation function of electron gas can be written as 

v(r) = < An(r 1 )An(r 2 ) > —£(7^—7*2) = 

2 , f ,. _ ,3 



yyexp (ipr /n)j£jh \ 2 , f = r 2 -r x . W 



< n > 1 J p ' K1 '■"'(2nn) 

Here the corner brackets stand for average, An(r) = n(f)— < n > is the depar- 
ture of electron density n{r) from its average value < n >, 

- Fermi distribution function, p- momentum, ep - Fermi energy. Note, that 
equation (1) is valid for any (not necessary isotropic) dispersion law e(p). 

Let us remind, that for e(p) = p /2m , equation (1) leads to the following 
expression for the CF[1]: 



3 n 



F 4 cos 2 ^,r > ft/ffF (3) 
2n pft^ n 
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(pF - Fermi momentum), which contains Friedel oscillations (FO) [2] with the 
wave number 2pp/h ■ 

One can easily determine the wave numbers of FO(ii the direction of radius 
vector f is given) also in the case of an arbitrary dispersion law. Indeed, the long- 
wavelength behaviour of the C-F(equation (1)) is determined by the singularities 
of the integral 

S(.P\\) = J n f d 2 p ± , P\\=pr/r, p±=p-p^/r. (4) 

as a function of p\\ . Obviously, S(p^) is nothing but the square of the section of 
the Fermi surface (FS) formed by a plane perpendicular to vector p and located 
at the distance p|| from the origin in p-space. Therefore the singularities we are 
interested in correspond to the tangencies of the FS to planes perpendicular to 
r. The nature of these singularities depends on the local geometry of the FS 
near tangency points. 

Therefore the wave numbers of FOAfaj are the distances in j?-space between 
planes tangent to the FS and perpendicular to r(see Fig.l). It is well-known, 
that the wave numbers Afcy correspond also to Migdal - Kohn singularities 
[3] in the spectrum of phonons propagating in the direction r/r, as well as to 
singularities in the longitudinal plasmon spectrum (cf. [4]) etcQ . Thus, the 
vectors Afcy r/r define the distinguished points in the reciprocal space; since the 
present paper is devoted to the spatial correlation of electrons, we shall study 
the effects that occur in the real space. 

One can also come to the same conclusions if one considers the Fourier 
component of the CF 

which is proportional to the volume of the intersection of the FS and its ana- 
logue, shifted by the vector hk. The singularities of this quantity (as a function 
of k) correspond to the tangencies of the FS to its shifted analogue, as well as 
to k = 0. We may conclude tentatively, (see section 2 below) that the way of 
damping of FO at r — > oo is determined by the local geometry of the FS at 
the tangency points. Thus, the behaviour of the CF, as well as many other 
electronic properties of metals (cf. [6]) is determined by the FS geometry. 

It turns out that one can describe the asymptotic behaviour of the CF at 
large r with the formula 

^ = -^4^E^iexp(^,r) . (6) 

ij 

2 Also, the quantities (Ak i ^)~ 1 are proportional to the periods of oscillations of the sound 
absorption coefficient in a magnetic field [5]. 
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Here r e =< n >~ 1 / 3 is the average distance between electrons, and indexes i,j 
enumerate the tangency points. Factors ^j( wmcn are > m general, functions of r 
as well as of r j r) depend on the local geometry of the FS at the corresponding 
tangency points and have a dimensionality of inverse length. Due to the central 
symmetry of the FS, expression (6) is real (obviously, Akij= —Akji). 

At i — j , Ak i j= and corresponding terms in the CF decrease monoto- 
neously as r increases. These terms are due to the singularity (namely, discon- 
tinuity of the first derivative) oiv{k) at k — 0. 

Usually, one can write down the validity condition for equation (6) as 

Akijr 3> 1 for any i ^ j . 

Note that the wave numbers Afc^ (as well as their total number) depend on the 
direction of radius vector r . 

So far we discussed the electron gas with an arbitrary dispersion law. The 
co-ordinate dependency of the CF of conduction electrons in a real metal is 
not exhausted by equation order to describe this case, one should take 

into account the influence of periodic potential of crystal lattice not only on 
the spectrum but also on the electron wave functions: wave functions that 
should be used now are the Bloch waves (instead of plane waves ). This topic 
is considered in detail in our previous article [7] . The qualitative results of this 
straightforward though somewhat cumbersome consideration are listed below. 

1. Due to the broken translation invariance, the co-ordinate dependency of 
the CFv(r 1 , r 2 ) is not reduced to the dependency only on the difference of its 
arguments r = r 2 — f\ . If the value of r is fixed, v{r\, r 2 ) is a periodic function 
of f\: 

v {r 1 ,f 2 ) = v{r 1 ,r + r 1 ) = v{r 1 +a,f '+ r^+a), (7) 

where a is any lattice period. 

2. Let us keep r= const and average v{f 1 , r 2 ) over f\. Denote the obtained 
quantity by v{r). Its asymptotic behaviour will be similar to the one described 
by equation (6). Wave numbers of FO are to be determined in the same way - 
the only difference is that now one should imagine the FS as a periodic surface in 
the reciprocal space (not restricted to the first Brillouin zone alone!). Therefore, 
the expression for v{r) should contain summation over reciprocal lattice vector 
6; terms with b ^ 0(Umklapp terms) oscillate with a wave number Akij + br/r. 

3. CF contains a term, whose oscillation wavelength is the inverse diameter 
of the FS. Note, however, that if the FS is open, there can exist such a region 
of radius vector r directions, where v(r) does not oscillate at all. 

Let us now consider the dependency of the CF on the local geometry of 
the FS. For the sake of simplicity, we shall restrict ourselves to the analysis of 
equations (1) and (6); this approach may be easily generalized for the case of 
conduction electrons in a real metal [7] . 
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2. The Effect of Local Geometry of Fermi Surface on Correlation 
Function. 



From equation (6), it follows that each term in the CF contains the factors Ai 
, which depend on the local geometry of the FS at each of the two contributing 
tangency points. The expressions for Ai(r) for various cases are listed below. 
We shall not describe here the calculations - one might find such a description 
in [7]. 

Let us keep the direction of radius vector r fixed and consider the factor 
, corresponding to the i th tangency point. 

(i) If it is an elliptic point, then 

A ell =2G- 1 ' 2 , (8) 

where G is the Gaussian curvature (product of the principal curvatures [8]) of 
theFS. 

(ii) For a hyperbolic tangency point, one obtains 

^hyp= ±2i I G T 1/2 > (9) 

where the sign depends on the direction of radius vector (whether r or — r). 

In these cases, A does not depend on r, and the amplitude of FO corre- 
sponding to a pair of elliptic and/or hyperbolic points of tangency, decreases at 
r — > oo as r~ 4 . 

(iii) The domains of elliptic and hyperbolic points on the FS are separated 
by the lines of parabolic points. Up to very few exceptions (Na,K,Cs,Rb 
and Bi), the FS of real metals are rather complex and do contain the lines of 
parabolic points [6,9]. 

In the generic case, one can choose the local coordinates £i , £2 on the FS in 
such a way, that the departure £3 of the FS from the plane, tangent at parabolic 
point, reads: 

« Cg-B& (10) 

(here we restrict ourselves to the leading order in x,y). We shall assume, for 
convenience, that B > and C > 0. 

If the tangency point is parabolic, then 

2IC( S inn/12)r(5/6) /6 /6 
^ipar- 3l /4 C1 /3 B i/2 n r ■ y 11 ) 

Here JC - elliptic integral, T - Euler's gamma-function; 

JC{simr/12)r(5/6) w 1.5. 

Thus we immediately notice, that if parabolic tangency points appear at a given 
direction of r then the amplitude of the corresponding term in the CF decreases 
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as r" 11 / 3 , i.e. slower than in the usual case. Due to the fact that parabolic 
points on the FS are not isolated, but form continuous lines, there exist cones 
of radius vector directions corresponding to the presence of parabolic tangency 
points^] . 

Let us now denote by A9 the angle between the normal at the parabolic 
point and the projection of radius vector onto the plane perpendicular to the 
line of parabolic points (Fig. 2). At AO < 0, there are no tangency points 
near the parabolic points. At A9 > 0, however, there exist both elliptic and 
hyperbolic tangency points and both of them approach parabolic point as A6 — > 
0. Therefore in the CF there appears an additional long- wavelength term which 
has the wave number 

g. S -3/2£7-l/2(^)3/2. 

Due to the vicinity of a parabolic point, equations (8) and (9) should be modified 
(one may not take into account terms quadratic in x, y only). Thus, for an 
elliptic point (Fig. 2, x < 0), we have 



5iC I 1 



16 ■ 

and for a hyperbolic point, we have 



A hyp a l B~ 1 /\3CA6)- 1/i {l+ - - — — 9 ( 1 (13) 



C 1 ' 2 5n 1 

2n3 1 / 2 (A6) 3/2 8 l) r> 

This distinction between elliptic and hyperbolic points is not so surprising. One 
has to remember that the singularities of v(k) corresponding to elliptic and 
hyperbolic tangency points are somewhat different. 
Both formulae (12) and (13) are valid at 

A8 » (r 2 /C) 1/3 (14) 

Suppose that at A9 < there exist 2M tangency points (here we are taking 
into account the symmetry of the FS) and, therefore, 

different wave numbers of FO. At A9= there appear two new (parabolic) 
tangency points, and the amount of FO wave numbers increases by Af + 2. 
Then at A6 > 0, A8<(r 2 / C) 1/3 the crossover occurs (the contribution of each 
parabolic tangency point turns into the sum of contributions of elliptic and 

3 If the FS contains parabolic points then the singularities of kinetic characteristics of the 
metal are enhanced for the corresponding directions of the quasimomentum vector [10,11]. 
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hyperbolic tangency points) and finally at AO S> (r 2 / C) 1 ^ 3 there are 2(M + 2) 
tangency points and (TV + 2)(N + 5) / '2 — 1 FO wave numbers. 

Let us consider the Gauss map of the FS(Fig. 3a). The Gauss map is a 
mapping of the surface onto the unit sphere, induced by the normal at each point 
of the surface (the direction of the normal gives a point of the sphere, thus giving 
the image of an original point of the surface) [8] . The number of tangency points 
changes by 4 as the point, representing the direction of the radius vector, crosses 
the line corresponding to the directions of the normals at parabolic points. This 
is the only way for this number to be changed. Suppose that for the direction 
of a radius vector lying within the area 1 (see Fig. 3a) the total number of 
the tangency points (on the whole FS) is 27V. Then for rjr situated on the 
lines BE or ED there are 2(J\f+l) tangency points (two parabolic points have 
been added). The number of tangency points for rjr lying in areas 2 and 4, on 
lines AE or EC, and in area 3 are given by 2{M + 2), 2{M + 3), and 2{N + 4 ) 
respectively. The angle 6 appearing on Fig. 2 is just the distance between point 
rj r and the image of the line of parabolic points on the Gauss map. 

(iv) Points at which the FS becomes flat [12] manifest another type of local 
geometry. The FS with such a point is a boundary case between convex surfaces 
and surfaces with a'crater'. Near this flattening point, the departure of the FS 
from the tangent plane is a fourth-order form of the local coordinates within 
the surface. If the tangency point is a point of such a kind, then 

A fl oc r 1/2 . (15) 

Let us assume for definiteness, that the departure of the FS from the plane 
tangent to the FS at the flattening point, is given in the leading order by 

= s£- {16) 

Then as the (elliptic) tangency point approaches the flattening point, the Gaus- 
sian curvature decreases and 

A Qll « 2 5 /6 5 -V2 £,-l/3 A 0-2/3 ^ ( J7 ) 

Here AO is now the angle between the radius vector and the normal to the FS 
at the flattening point. Formula (17) is valid at 

AO > D 1 /4 r -3/4 ^ 

where the intermediate region is the region of crossover between the dependen- 
cies (15) and (17) (see Fig. 36). 

(v) Suppose that the lines of parabolic points on the FS cross at some point. 
In real metals such a crossing occurs rather often; for example, the FS of the 
form shown in Fig. 4 is encountered in Mo, W , paramagnetic Cr[6, 9, 11]. 

The departure of the FS from the plane tangent to the FS in this point of 
crossing, is given by a third-order form and we have 



G 



Aqi oc r 1/3 



(19) 



In fact, this is a flattening point of another kind. The important feature of 
flattening tangency points of any type is that they correspond to isolated points 
on the Gauss map (see Fig. 3). 

(vi) Let us assume that the FS contains finite flat elements. If the radius 
vector is parallel to the normal to such an element of the FS, it leads to the 
contribution of this "tangency area" in the CF: 

^plane x ^plane r ' (^) 
where <Spj ane is the square of this flat element. 

In fact, the FS of such a kind can hardly be stable: in this case, Peierls 
transition [13] with the transformation of the lattice structure should become 
energetically favourable (at least at sufficiently low temperatures). It is easy to 
notice that there should be some correlation between extremely slow damping 
of FO as r — > oo(here v oc r~ 2 cos (Z\fcp^ ane r)) and the possibility of Peierls 
transition. 

At this point we finish our consideration of the local geometry of the FS. Note, 
that we have arrived to the remarkable conclusion: the exponent, which de- 
termines the damping of FO as r — > oo ; depends on the local geometry at the 
points of tangency. Thus, one might build up a hierarchy of types of the FS 
local geometry with respect to the long-distance behaviour of the CF: 

1. Spherical FS: Amplitude of FO decreases as r -4 ; this exponent is given, 
generally, by elliptic or hyperbolic tangency points and thus corresponds to a 
generic direction in aSD metal. 

2. Cylindrical FS: Amplitude of FO decreases as r -3 . This is the case of 
a2D metal, or the case of a toroidal FS cavity (see sect. 5 of Ref. [7])Q . 

3. Flat FS. Amplitude of FO decreases as r -2 (one could imagine it as aFS 
of an ID metal). 

From our consideration, it follows that parabolic points (along with crossing 
points of the lines of parabolic points) should be placed "between" spherical 
and cylindrical FS. 

3. Some Applications. 

Let us consider a metal near the electron topological transition (ETT) point 
[14]. The ETT is the restructuring of the FS which occurs when ep= e cr (e cr 
- the value of electron energy that corresponds to Van Hove singularity). De- 
pending on whether this singularity is due to an extremum or to a saddle point 
of the electron dispersion law, a new cavity of the FS appears (disappears), 
or a "neck" connecting two FS cavities is formed (disrupted). The difference 

4 Point at which the FS becomes flat also gives the r~ 3 law (see Eqn (15)). 
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z = 6f— £cr depends on the applied pressure (or, for example, impurity concen- 
tration), so that the ETT can be realized experimentally (see Ref. [15] and 
references there). The ETT causes distinctive singularities of thermodynamical 
and kinetic quantities [14-16]. 

Near the ETT point, long- wavelength terms appear in the CF due to the 
existence of small diameters of the FS. If the ETT results in the appearance of 
a new cavity of the FS at z > , then 

"ETT(*0 * -<^?-#- ^ff^^ cos 2 {(^Eti miCOB'Of/nUSl) 

Here mi, m 2 , m 3 are effective masses and the axes x\, £ 2 , f 3 (where Xi— r cosOi) 
are the main axes of the tensor of inverse effective masses, so that the electron 
dispersion law near the point of extremum can be written as 

2m 1 2m 2 2m 3 

If (instead of an extremum) one has a saddle point of the dispersion law, one 
should write instead of (22): 

(we assume that the FS near the extremal point Ap = possesses the rotational 
symmetry with respect to the p 3 axis and is also symmetric with respect to the 
AP3= plane). In this case, the "neck" of the FS, which exists at z > 0, is 
ruptured at z < 0(Fig.5). Therefore, at | z |<C e F there exist long-wavelength 
(Ak<(2m | z \) 1/2 /h < k F ) terms in the CF. The issue of significance here is 
that the angular dependence of these terms "rotates" by ir/2 as the sign of z 
changes. 

Indeed, for example, at z < and r [[ p 3 , due to the small distance between 
the two sheets of the FS, there exist long-wavelength component of the FO 

1/2 

with corresponding Ak — {2m ^ \ z \) /% and amplitude 

m\ | z | r i /{2U 2 7r 4: < n > mn). 

There are no long- wavelength oscillations at f _L p 3 . Vice versa, for z > , the 
long- wavelength oscillations are absent at r || p 3 . For z > and r _L p 3 , due 
to small diameter of the FS "neck", the following long- wavelength term in the 
CF appears: 

"ETTW ~ - 4 ^ m "^ 4*2 sin2 ( r {2m ± z) 1/2 ) at r ± p 3 . 

At (3z > , the fourth-order term in {23) gives rise to the appearance of the 
lines of parabolic points on the FS near p — p cr [16]|^] (see Fig. 5). Thus for 

5 Note that parabolic points on the FS themselves are not peculiar from the point of view 
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some directions of r forming angle 9pai with ps axis, the amplitude of the 
long-wavelength term in the CF behaves as r - n / 3 — r -4 r 2/3 ms ^ ea( j f r -4_ j n 
the neighbourhood of these peculiar directions, there exists an additional small 
wave number of FO, which is due to neighbouring tangcncy points pairs (elliptic 
and hyperbolic). At 9= 6*par , each pair "sticks together", giving the parabolic 
tangency point and then disappearing. 

On the whole, the ETT, being a local "event" in reciprocal space, in co- 
ordinate space results in dramatic restructuring of the long-wavelength terms 
in the CF. 

Obviously, the knowledge of the CF is useful when one considers various 
collective phenomena in metals: screening, structure of spin glasses etc. Let 
us mention for example the RKKY interaction [17-19]. This indirect exchange 
interaction between the site (or impurity) spins is caused by exchange interaction 
between the localized spins and conduction electrons: 

H=(J/ <n>)Y / ^S l S(r-r l ). (24) 

i 

Here J - exchange constant, a - Pauli matrices, Si - spin of the ion located in 

the site n. The RKKY exchange integral reads 

jRKKYir) = ~2{^f ! n ^^e^[i{q-q')r/n^f,{25) 

and its Fourier component, 

JfcKKY(*) = -^)'J ^(& - W 

This kind of integral appears also in studying of the electron-phonon interaction. 
It has singularities at the same points as does the CF(except k = 0). These 
singularities are called Migdal - Kohn singularities [3] if the electron velocity 
vectors in the two contributing tangency points are antiparallel. On the other 
hand, if they are parallel, then the singularity is called Taylor singularity [20] 
and these two cases should be treated separately. 

It turns out, that at least for the great majority of cases there exists a quite 
simple relationship between the asymptotics of v(r) and ^RKKy(^)' Indeed, 
for a fixed direction of the radius vector 

■ / mcKY ocri/(Afc) ( r )' Ak *° W 

where J^^y anc ^ ^ Ak H r ) stand for the terms in Jp^xKY anc ^ the respec- 
tively, oscillating with the wave number Ak. 

It is well-known, that magnetic order induced by the RKKY interaction 
(namely, the helicity vector) depends on the geometry of the F5[18]. Thus, the 
helicity vector is expected to be determined by the small diameter of the F5[21]. 

of electron dispersion law. Generally, the lines of parabolic points may have nothing to do 
with ETT or Van Hove singularities (see sect. 2 above). 
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Therefore the knowledge of the FS and the CF seems to allow one to make some 
predictions about the magnetic arrangement. We expect, in particular, that the 
ETT may manifest itself in changes of magnetic ordering. 

4. Conclusion. 

In this paper, we neglected temperature effects, electron scattering and 
electron-electron interaction. The first two effects lead to the appearance of 
an exponentially decreasing factor in the CF as r — > oo. The third (i.e. Fermi- 
liquid interaction), probably leads to the renormalization of the coefficients (see, 
for example, [22]). 

Anyway, the dependence of the asymptotic behaviour of the correlation func- 
tion on the FS geometry (including the local geometry) survives after taking 
these effects into account. This dependence should affect collective phenomena 
in metals. Indeed, one has to know the CF to construct theories of alloying, 
screening, exchange magnetism etc. We have investigated some features of the 
CF that seem to be interesting from this point of view. Among them are: 

1. The existence of the cones of peculiar directions, corresponding to the 
lines of parabolic points on the FS. For these directions of the radius- vector, the 
CF shows some specific features: 

(i) the attenuation of FO is determined by a factor of r~ n / 3 instead of usual 
r~ 4 ; we doubt if this small difference has any observable effect. 

(ii) the number of FO periods changes as the direction of the radius vector 
crosses these cones. 

(iii) On one side of such a cone, one of FO periods appears to become large 
(it goes to infinity as the direction of r approaches the cone) . 

2. The existence of the isolated peculiar directions. These correspond to the 
crossings of the lines of parabolic points on the FS* (flattening points). 

3. The restructuring of the CF due to electron topological transition. 

We have also provided one example showing how an asymptote of the RKKY 
exchange integral depends on the asymptotic behaviour of the CF. 

It is a pleasure to thank M.Yu.Kagan and L.P.Pitaevsky for fruitful dis- 
cussions. We are also grateful to A. Prakash for assistance in preparing this 
manuscript for publication. 
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Figure Captions 

Fig.l Rather simple FS(" dumb-bell" ) for the radius vector r \\f\ generates 8 
different wave numbers of Friedel oscillations (due to the symmetry, 

Afci2= Ak$6, A&i3= Ak i6 , Ak u = Ak 2 5= Ak 36 , 
Afci 5 = Ak 2 Q, Afc 2 3= ^&45and Afc 2 4= ^35)- 

Circled numbers enumerate the tangency points for this case. The total number 
of wave numbers depends on the direction of r: for r \\ r 2 there is the only wave 
number Ak. 

Fig. 2. Parabolic point on the FS^placcd at the origin of local co-ordinates 
^17^27^3 7 where £3 is parallel to the normal). For a small value of angle 9 > 0, 
there are two tangency points in the vicinity of the parabolic point: elliptic 
(with £1 < 0) and hyperbolic (with £1 > 0), and, correspondingly, the small 



11 



wave number Ak oc 9 3 ^ 2 . As 6 — > +0, these tangency points approach each 
other; at 0= they coincide at the parabolic point; at 6 < 0, there are no 
tangency points in the vicinity of this parabolic point (see also Fig. 3a). 

Fig. 3. Gauss maps of the elements of the FS. 

3a. The image of the part of the FS with two lines of parabolic points AC 
and BD intersecting in the point E. Each point in areas 2 and 4 is the image of 
the two points from the original part of the FS, each point in the area 3 - the 
image of four points and each point in area 1 does not have any original there. 
Shaded areas correspond to crossover (the inequality (14) is broken there). This 
obviously is only one of the two kinds of intersections of the lines of parabolic 
points. 

3b. Gauss map of the part of the FS containing the point at whichthe FS 
becomes flat. The number of tangency points is constant for r/r laying in this 
region of directions. The shaded area corresponds to crossover (inequality (18) 
is broken there). 

Fig. 4 Crossing of the lines of parabolic points on the 

4a. Fragment of the FS(" octahedron") in a metal belonging to the molyb- 
denum group. 

4b. Crossing of lines of parabolic points on this cavity of the FS (view from 
above) . 

Fig. 5 Disappearance (a) and appearance (b) of the lines of parabolic points 
(marked by the thick lines) during the "neck" formation [16]. 
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